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Motivated by recent huge interests on graphene sheet and graphite as "relativistic" systems, a 
BCS superconductivity in a quasi-(2+l)-dimensional relativistic model is investigated. The intra- 
layer particle dynamics is described by a (2+l)-dimensional Gross-Neveu-type four-body contact 
interaction model, while inter-layer particle motions will be caused by a hopping term with a transfer 
parameter t. Especially, we examine the effects of non-vanishing t, chemical potential /x, and a mass 
parameter m which will give a gap at a conical intersection point of the relativistic band dispersion, 
in the BCS s-wave ( scalar ) superconducting gap function A. 
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The honeycomb lattice structure has been obtained remarkable interests recently in condensed matter physics. 
For example, graphite and a graphene sheet have the honeycomb lattice symmetry, and they show various interesting 
physical properties [1-22]. It is well-known fact that, a band structure of a system of honeycomb symmetry has conical 
intersections between the valence and conduction bands, and the band dispersions can be described by a Dirac model. 
Graphite can be regarded as a system where (2+l)-dimensional layers are stacked in a specific direction. Toward some 
efficient design for electronic substances, it is desirable for us to find simple but powerful tendencies/laws ( especially 
as a function of inter-layer transfer t, which is approximatedly given by an overlapping integral between orbitals of 
two layers ). A study by employing a simple field-theoretical model seems suitable for us for this purpose. 

Toward an examination on superconductivity in the honeycomb lattice symmetry, the author studied a BCS-type 
theory in a single layer system ( the exact (2+l)-dimensional case ) in Ref. [7] as a first step. Recently, Ref. [5] 
investigated a kind of BCS-type superconductivity in a lattice model of graphene sheet. It is interesting for us to 
investigate some effects coming from an inter-layer interaction ( particle hopping between two nearest neighbor layers 
) in relativistic BCS superconductivity. Of course, such motions of particles by hopping between layers surely exist 
in real substances, except the case of an ideal mono-layer situation. Therefore, we will investigate a BCS theory in a 
quasi-(2+l)-dirnensional system where x and y directions are homogeneous, while particle motions of z-direction are 
described by a Hubbard-like hopping ( lattice ) model in this paper. 

The energy spectra of particles near a conical intersection point with having a hopping term with a transfer 
parameter t perpendicular to two-dimensional layers will be given as follows: 



E { l\k) = ±\kx\ -<cosfc 3 , k x = (ki > k 2 ). 



(1) 



We will seek a model Lagrangian which will give the energy spectra of (1) from its kinetic term, while it also has a 
BCS-type interaction between particles. The candidate for us will be found in the following form: 
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( a similar model was discussed in Ref. [23] ). Ck has both the usual Dirac kinetic term in (2+l)-dimensions and the 
inter-layer hopping terms. The abbreviations for coordinates and derivatives inside the layer ( i.e., perpendicular to 
the third, stacking direction ) are defined by 



x± = (xi,x 2 ), d J _ = (d 1 ,d 2 ), 
while, we choose the following definition of the three-demensional gamma matrices: 



(5) 



7° = cr 3 , 7 1 =icn, 7 2 = ier 2 , 7 ± = (7 1 ,7 2 ). 



(G) 
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ip is a two-component spinor. The Fourier transform of the fermion field in our model becomes 

4>(x ,x ± ,n) = J iP(ko,k ± ,h)e-* kaX °+^ x ^ k3n , (7) 

/<••■>- /$/&/>■•>• 

The integration of k 3 is restricted to the region — tt < k 3 < tt, namely the first Brillouin zone. Thus, the Fourier 
transformed expression for the kinetic Lagrangian is found to be 

Ck — rp(j°{ko + n + t cos/c 3 } + j 1 - ■ kj_ — m — ucos k^j ip. (9) 

The diagonalization of Ck at /x = m = u = gives (1) as its eigenvalues. Our Lagrangian C becomes a (2+1)- 
dimensional Gross-Neveu model [7,24,25] at the zero transfer limit t — > 0, u — > 0. ucosfc 3 is a mass with an 
alternation ( modulation ) in the fc 3 -direction. 

In this paper, we only consider scalar-type fermion pair condensates a ( fermion-antifermion ) and A ( fermion- 
fermion ). From possible functional ( mathematical ) structure of gap equations, this choice is relevant for considering 
BCS superconductivity [7]. We introduce the method of auxiliary fields of composites to our Lagrangian: 

C = C K - (rrpip + ^4>C^ T + ^-ip T C- V ^ 



2 r r 2 r T 2G 2G 
= -±(a* + \A\*)+hM*. (10) 

Here, C is the charge conjugation matrix, and we have taken the four-component Nambu notation [26] of the following 
definition: 

•<*>-( $2> )• •«-(**>.*->». <"> 

where, 

Mn = J°{ida +/i + tcosfc 3 } + 17 • d±_ — a — m — u cos k$ , 
X 12 = A(x)C, MjisSfiJC- 1 , 

A^22 = 7° T {i<9o — /U — tcosfc 3 } + i(7 ± ) T • <9± + ct + m + ucosfc 3 . (12) 
Hence, the partition function Z considered in this paper is defined as follows: 

Z = J V^VtpVaVAVAexpU^ / dx dx 1 dx 2 C . (13) 

The effective potential becomes 

cr 2 |A| 2 /■ 

where the determinant is evaluated to be 

tr f \ndctM = [ ln[fc - E+] [k + E+] [k - E_] [k + E_] . (15) 
The eigenvalues appeared in V e /f will be obtained in the following forms: 



E±(k) = y (^fc 2 + [a + m + ucosfc 3 ] 2 + [^ + icosfc 3 ]) + |A| 2 . (16) 

We employ the finite-temperature Matsubara formalism for our convenience of /co-integration [27-29]. The zero- 
temperature theory is converted into the Matsubara formalism by the following substitution in our theory: 



/^^X!^' k o^iu n , P=7f, (T; temperature), ui n = ( 2n + 1 ) 7r ? (n = 0,±l,±2,- 



(17) 
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After performing the discrete frequency summation, we get 

a 2 |A| 2 <- K ^ 2 



eff 



1 0/3(^—1^—)'— ))■ 



2G 2G 

where, A is a two-dimensional momentum cutoff: 

The gap equations, namely the stationary condition of V e ff will be obtained as follows: 

dVe/f = o_ _ r A fk^ r dh 

da G J (2n) 2 J_ 7T 2n 

and 

dV eff _ A f A d 2 k ± f n dk 3 



— t tanh + -r— tanh -£L 
9cr 2 + da 2 



= 
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Here, several derivatives in the gap equations become 

dE± \Jk\ + [a + m -\- u cos k 3 ] 2 =p [/z + i cos k 3 ] 
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Our gap equations have model parameters as G ( coupling constant ), A ( two-dimensional momentum cutoff ), 
m ( gap parameter at the conical intersection ), u ( a mass parameter alternating to the inter-layer direction ), t 
( transfer matrix between two nearest neighbor layers ) and /i ( a chemical potential ). We choose the unit A = 1 
throughout this paper. The coupling constant G will be chosen to satisfy the condition a, A -C A. In this paper, we 
do not consider a coexistent phase of a ^ and A^O, and study phases of both a ^ 0, A = and a = 0, A ^ 
separately. 

First, we investigate the critical coupling for dynamical generation of a. At m — u — A = T = 0, the self-consistent 
gap equation for a becomes like 

and this equation gives the critical coupling G cr as 

G C r = ■ (25) 

Hence, one cannot find any effect of t in G cr atm = !i = A = T = 0, namely the inter-layer hopping effect of t ( 
and also fi ) disappears in G cr at this case. If we consider the case A ^ with m = u = T = 0, G cr will obtain an 
expression which includes an effect of t, though we have to solve the gap equations of a and A simultaneously for our 
confirmation of coexistence of dynamically generated a and A. This is an interesting problem, though it is beyond 
scope of this paper. 

Numerical results from the gap equation for A are given in Figs. 1, 2 and 3. Figure 1 shows the BCS gap function 
A at a = as a function of the inter-layer transfer parameter t. We set parameters as A = 1 (unit), G = 4, 
T = u = fj, = 0(fj, = corresponds to the "Dirac" vacuum state ). The four examples of m = 0,0.1,0.2,0.3 arc 
depicted in this figure. When t ( m ) increases, A increases ( decreases ). There is no non-trivial solution for A at 
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the cases of t = 0, m = 0.3 and t — 0.1, m — 0.3, and this fact indicates the existence of a "critical transfer" ( and 
also, a "critical coupling" ) for dynamical generation of finite A in our theory. 

Figure 2 gives the BCS gap function A at a — as a function of t. The parameters A, G, T and u are set as same 
as Fig. 1, while we choose \x = 0.2. Increasing of fi causes the enhancement of A. 

Figure 3 shows A at a — shown as a function of t. The parameters A, G, T and u are set as same as Figs. 1 
and 2, and chemical potential is chosen as \x — 0.4. In this figure, A decreases when t increases. Hence, i-dependence 
of A is not trivial, and it is affected by [i. It is also the case that, when m is small enough, the difference of choices 
of numerical values of fj, disappears at large t ( effect of fj, is washed out by a large t ) . See the differences of energy 
scales of Figs. 1, 2 and 3. 

Temperatre dependence of gap function A is shown in Figs. 4 and 5. We choose three values of /i in Fig. 4, while 
the three examples of t are chosen in Fig. 5. In both of these figures, the ratio A(T = 0)/T c ( T c ; critical temperature 
) becomes 1.7 ±0.1, i.e., the BCS universal constant A(T = 0)/T c = 1.76 is almost satisfied, though it seems the case 
that there are tiny but finite deviations from the constant. The gap function always show the character of second 
order phase transition. 

We summarize our conclusions as the following tendencies: (1) When m increases, A decreases. (2) When n 
increases, A increases. (3) A depends on t non-trivially, and its dependence on t is qualitatively changed by choices 
of numerical values of \i. 

If we consider a real substance like graphite, fi corresponds to a charge doping concentration to the system. Our 
results indicate that a charge doping and an inter-layer coupling do not always give cooperative effects to an amplitude 
of a BCS gap function. The results of this paper can also be used to consider superconductivity of MgB 2 [30] which 
also has the honeycomb lattice symmetry [7]. It is very interesting for us to use parameters obtained by an ab initio 
band structure calculation or a molecular orbital calculation for getting numerical values of t, so forth. 



* Present Address; Department of Physics, University of Texas at Austin,tadafumi@physics. utexas.edu 

1 G. W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

2 F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

3 J. Gonzalez, F. Guinea and M. A. H. Vozmediano, Phys. Rev. Lett. 77, 3589 (1996), Phys. Rev. B63, 134421 (2001), J. 
Gonzalez, Phys. Rev. Lett. 88, 076404 (2002), N. M. R. Peres, A. H. Castro Neto and F. Guinea, Phys. Rev. B73, 195411 
(2006), A. H. Castro Neto, F. Guinea and N. M. R. Peres, Phys. Rev. B73, 205408 (2007), E. V. Castro, K. S. Novoselov, 
S. V. Morozov, N. M. R. Peres, J. M. B. Lopes dos Santos, J. Nilsson, F. Guinea, A. K. Geim and A. H. Castro Neto, Phys. 
Rev. Lett. 99, 216802 (2007), J. Nilsson, A. H. Castro Neto, F. Guinea, and N. M. R. Peres, arXiv:0712.3259, F. Guinea, 
B. Horovitz and P. Le Doussal, Phys. Rev. B77, 205421 (2008), F. Guinea, arXiv:0805.3908. 

4 A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov and A. K. Geim, arXiv:0709.1163. 

5 B. Uchoa and A. H. Castro Neto, Phys. Rev. Lett. 98, 146801 (2007). 

6 D. V. Khveshchenko, Phys. Rev. Lett. 87, 246802 (2001), Phys. Rev. B74, 161402(R) (2006), Phys. Rev. Lett. 97, 036802 
(2006), Phys.Rev. B75, 153405 (2007). 

7 T. Ohsaku, Int. J. Mod. Phys. B18, 1771 (2004). 

8 V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 146801 (2005). V. P. Gusynin, V. A. Miransky, S. G. Sharapov and 
I. A. Shovkovy, cond-mat/0612488, arXiv:0801.0708, arXiv:0806.2136. E. V. Gorbar, V. P. Gusynin, V. A. Miransky and I. 
A. Shovkovy, arXiv:0806.0846. E. V. Gorbar, V. P. Gusynin and V. A. Miransky, arXiv:0710.3527. 

9 C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005), C. L. Kane, Int. J. Mod. Phys. B21, 1155 (2007). 

10 N. A. Sinitsyn, J. E. Hill, H. Min, J. Sinova and A. H. MacDonald, Phys. Rev. Lett. 97, 106804 (2007). 

11 K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 96, 256602 (2006), Phys. Rev. Lett. 98, 076602 (2007). 

12 K. Yang, S. Das Sarma and A. H. MacDonald, Phys. Rev. B74, 075423 (2006). 

13 Y. Barlas, T. Pereg-Barnea, M. Polini, R. Asgari and A. H. MacDonald, Phys. Rev. Lett. 98, 236601 (2007). 

14 H. Min, G. Borghi, M. Polini and A. H. MacDonald, Phys. Rev. B77, 041407(R) (2008). 

15 M. Polini, R. Asgari, Y. Barlas, T. Pereg-Barnea and A. H. MacDonald, Solid State Commun. 143, 58 (2007). 

16 M. Polini, A. Tomadin, R. Asgari and A. H. MacDonald, arXiv:0803.4150. 

17 H. Min and A. H. MacDonald, Phys. Rev. B77, 155416 (2008), arXiv:0806.2792. 

18 D. Xiao, W. Yao and Q. Niu, Phys. Rev. Lett. 99, 236809 (2007). 

19 Y. Hatsugai, T. Fukui and H. Aoki, Phys. Rev. B74, 205414 (2006). 

20 M. Kohmoto, arXiv:0709.2213 

21 P. Ghaemi and F. Wilczek, arXiv:0709.2626. 

22 D. L. Bergmann and K. Le Hur, arXiv:0806.0379. 

23 V. M. Loktev, R. M. Quick and S. G. Sharapov, Physica C314, 233 (1999). 

24 K. G. Klimenko, Z. Phys. C37, 457 (1988), C54, 323 (1992). 



■5 



FIG. 1: A shown as a function of t. We set parameters as A = 1 (unit), G = 4, T = u = a — n = ( n — corresponds to the 
"Dirac" vacuum state ). The four examples of m = 0, 0.1, 0.2, 0.3 are depicted. 

FIG. 2: A shown as a function of t. We set parameters as A = 1 (unit), G = 4, T = u = a — and fi = 0.2. The four examples 
of m = 0, 0.1, 0.2 are depicted in this figure. 

FIG. 3: A shown as a function of t. We set parameters as A = 1 (unit), G = 4, T = u = a — and \i = 0.4. The four examples 
of m — 0, 0.1, 0.2, 0.3 are depicted in this figure. 

FIG. 4: A as a function of temperature T. We have set parameters as A = 1 ( unit ), G = 4, u — a — 0, m = 0.1 and t = 0.2. 
Three examples of n — 0, 0.2, 0.4 are shown in this figure. 

FIG. 5: A as a function of temperature T. We have set parameters as A = 1 ( unit ), G = 4, u = a = 0, m = 0.1 and fi — 0.2. 
Three examples of t = 0, 0.2, 0.4 are shown in this figure. 
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